# Lesen der Gebäudedaten
gwr <- read.csv(paste0("https://www.web.statistik.zh.ch/ogd/daten/",
"ressourcen/KTZH_00002022_00004064.csv")
) |>
dplyr::rename_all(tolower)
# Umwandeln in ein sf-Objekt -> Geodatenobjekt
gwr_sf <- sf::st_as_sf(
gwr,
coords = c("e.gebaeudekoordinate", "n.gebaeudekoordinate"),
crs = 2056
)
# interaktives Darstellen
mapview::mapview(head(gwr_sf, 100))
gemeinden <- sf::st_read(
paste0("https://maps.zh.ch/wfs/OGDZHWFS?SERVICE=WFS&VERSION=2.0.0&",
"Request=getfeature&TYPENAME=ms:ogd-0095_arv_basis_up_gemeinden_seen_f&",
"outputformat=geojson")
) |>
dplyr::rename_all(tolower)